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Abstract 



An investigation of the kinematics of a rotating relativistic plasma stream 
in the perpendicular rotator model of the pulsar magnetosphere is presented. 
It is assumed that the plasma (ejected from the pulsar) moves along the pulsar 
magnetic field hues and also corotates with them. The field lines are considered 
to be radial straight lines, located in the plane which is perpendicular to the 
pulsar rotation axis. The necessity of taking particle inertia into account is 
discussed. It is argued that the " massless" (" force- free" ) approximation cannot 
be used for the description of this problem. The frame selection is discussed and 
it is shown that it is convenient to discuss the problem in the noninertial frame 
of ZAMOs (Zero Angular Momentum Observers). The equation of motion and 
the exact set of equations describing the behaviour of a relativistic plasma 
stream in the pulsar magnetosphere is obtained. The possible relevance of 
this investigation for the understanding of the formation process of a pulsar 
magnetosphere is discussed. 
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Introduction 



In connection with several astrophysical problems, e.g. the theoretical de- 
scription of pulsar magnetospheres, the question of investigation of plasma 
ejection from a rotating source (e.g. a star) along its magnetic field lines is 
arising. It is most important to investigate the behaviour of this stream near 
the pulsar surface. 

While investigating the pulsar wind, the change of particle momentum is 
usually neglected. The equations are then formulated in the "massless" 
approximation m—0 (Mestel, 1973; Endean, 1974), in which the particles are 
moving along the pulsar magnetic field lines like photons moving along geodetic 
lines. It seems to us indubitable, that the inertia of the particles, disregarded 
in a large number of previous papers, must have a great influence on the pulsar 
magnetosphere formation process. The importance of taking into consideration 
the inertia of the particles was understood a long time ago and treated as a 
perturbation to "massless" approximation (Scharlemann, 1974; Henriksen and 
Norton, 1975). It seems to us that giving up the "massless" approximation 
not only shows the essence of the problem more adequately, but also simplifies 
the problem from the mathematical point of view. 

We use the perpendicular rotator model of the pulsar magnetosphere and 
treat only the polar cap. It is assumed that the pulsar magnetic field lines are 
the radial straight lines. In our paper is discussed the case, when the plasma 
particles corotate with the pulsar magnetic field lines. We also assume that 
the E\\ — 0, i.e. the electric field component along the magnetic field lines is 
zero. 

We will choose the simplified geometry, discuss the frame selection and the 
derivation of the equation of motion of a relativistic stream. We will argue that 
the "massless" ("force-free") approximation cannot be used for the description 
of this problem and discuss the difference of our statement of the problem from 
those discussed before by a number of authors (Mestel, 1973; Endean, 1974; 
Henriksen and Norton, 1975; Michel, 1969; Kennel et al., 1983 and others). 

We will obtain the equation of motion of a relativistic plasma stream and 
the exact set of equations for the description of the behaviour of the stream. 

Basic Considerations 

Let us consider the simplified geometrical model, when the pulsar rotation 
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axis is perpendicular to its magnetic momentum. The magnetic field lines are 
considered as radial straight lines, located in the plane which is perpendicular 
to the pulsar rotation axis. If we discuss the processes in the region which 
has the linear size d<^Rc, where Rc is the curvature radius of the magnetic 
field lines, this approximation is justified. The ejected particles not only move 
along the radius r, but also coroatate with the pulsar magnetosphere because 
the field lines are corotating with the plasma: 

E+[VB]=0. (1) 
Here and below we use geometric units c=G=l. 

The condition (1) means, that because of the rotation of the magnetic field 
B, a field E is generated, which forces plasma particles to corotate with the 
magnetic field lines with linear velocity: 

(2) 

— * 

where Q is the angular velocity vector of the star. 

In the "massless" approximation, in an inertial frame there are no other 

— * — * — * 

forces in the equation of motion than the Lorentz force e{E + [VB]). The 
centrifugal force is neglected. However, by condition (1) the Lorentz force 
vanishes and ^=0. According to Rylov (1989) the equation ^=0 in approxi- 
mation m—0 can be justified in two cases: p^O and 7— >oo, or finite 7 and 
p—0. The first case is called by Rylov "force- free" approximation and the se- 
cond one — "massless" approximation. But to disregard the centrifugal force 
it must be shown under which conditions this is a valid approximation. Below 
we will show, that disregarding the centrifugal force is equivalent to giving up 
corotation. In order to understand the physical reason for corotation, let us 
consider a well known example: let us make a conductor move across a mag- 
netic field with a constant velocity V. Then, the Lorentz force e[yS] will act 
in different directions on the charges of different signs of the conductor which 
will cause a charge separation, i.e. generation of the electrostatic field. Due 
to the separation of charges, some energy of the source making the conductor 
move across a magnetic field lines with constant velocity (this velocity satisfies 
the condition (1)), must be expended. An analogous situation holds in the 
case when an observer moves with respect to the conductor with the velocity 
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V. The condition (1) allows also such an interpretation: the motion of the 
particles (of any sign) perpendicularly to the magnetic field is possible if an 
electric field E (EJ^B) exists and the velocity of this motion is equal to the 
electric drift velocity V=^^^. 

Now let us assume that the outer source makes the conductor move with an 
acceleration. Then it is evident, that the work of the source is expended not 
only on the charge separation, but also on the acceleration of the conductor. 
Therefore, except of the Lorentz force, there will be an additional force F in 
the equation of motion, but due to the condition (1) it will take the following 
form: 

dt 

Returning to our problem we conclude: the pulsar rotation makes the mag- 
netic field fines rotate, this rotation generates the electric field, which is ori- 
ented perpendicularly to the magnetic field lines and because of the condition 
(1) makes the charged particles (ejected from the star) corotate with the mag- 
netic field lines. During this motion, the particle velocity changes its direction, 
therefore a force is acting on the particles. This force in our case is directed 
along the magnetic field lines and is not compensated. 

Considering analogous problems in classical mechanics, for example the 
motion of a bead located in a rotating tube, the reaction force appears in the 
equation of motion. This force makes the bead not only rotate together with 
the tube, but also accelerates it along the tube. 

In our problem, according to the condition (1) the magnetic field lines and 
the generated electric field act like the rotating tube for the charged particles. 
It is more convenient to find the form of the reaction force in a noninertial 
frame, which rotates together with the particles. The situation described above 
is very similar to the problems connected with the description of particle mo- 
tion in curved space-time, in particular in the gravitational field of the black 
hole — in the Schwarzschild or Kerr metric, because according to the Einstein 
principle of equivalence, locally we cannot distinguish gravitation from non- 
inertiality. Using the principle of equivalence, the (3+1) formalism is suitable 
for our problem. The (3+1) formalism is described in Thorne et al. (1988). In 
this formalism instead of four- vectors, three spatial and one time coordinates 



5 



are used. 

The problem of the motion of particles in the pulsar magnetosphere can 
be solved in the local-inertial frame of the observers, who are measuring the 
physical quantities in the immediate vicinity of themselves. They are called 
the Zero Angular Momentum Observers (ZAMOs). ZAMOs cover the whole 
space-time. Each of the observers uses its proper time r, which is different from 
the universal time t. We must mention that the proper time of the observer 
who moves together with the particle, is different from the proper time of 
ZAMOs r. 

The transition from the inertial frame to the frame, which is connected with 
the pulsar magnetosphere rotation can be done by following relation: 

t=t' , (p=Qt , r=r' , z=0 . (3) 
The metric, describing such a frame has the following form: 

dS^=-{l-n\^)df+dr\ (4) 

^Prom (4) it follows, that the interval of the proper time of ZAMOs is 
connected with the universal time interval as 



dr—adt, (5) 

where a—y/l — Q^r^ and it is called the "lapse function". 

A photon, propagating from the periphery to the center of rotation, from 
observer to observer would be "reddening". The reason is that the rate of the 
time lapse for each of the ZAMOs differs from those for the others and also 
differs from the lapse of the universal (laboratory) time. The lapse function a 
not only connects the proper time of ZAMOs r with the universal time t, but 
also expresses the gravitational potential, i.e. the inertial force. In the frame 
of ZAMOs the particle has the gravitational acceleration 

9— — . (6) 
The equation of motion of the particle in the frame of ZAMOs takes the 
following form: 

f=^g+^{E+[VB]), (7) 
dr m 
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where 'y=—^=^ is the Lorentz-factor and V=-^^ is the velocity of the particle, 
pis the dimensionless momentum, p— > — . 

We must mention, that the equation of motion contains the acceleration 
g from (6) only in the frame, for which d(/9'=0. In the generalized rotating 
frame, where d^p'^O, according to the (3+1) formalism, the lapse function a=l, 
but instead of the acceleration (6) the "gravimagnetic" acceleration appears 
(Thorne et al. 1988), which must be written in the equation of motion. This 
equation differs from (7) only by the rcdcsignation of the angular velocities 
and, of course, describes the same process and gives the same results as (7). 

While considering the hydrodynamical problem, according to the (3+1) 
formalism, instead of the equation of motion for a single particle, the energy- 
momentum tensor conservation (the law of force balance) is used. Also, using 
the particle consertvation law (the continuity equation) and neglecting the 
hydrodynamical stream pressure, the equation, which describes the stream 
motion takes the following form: 

i|+(VV)p=-.^+^(i+(^i5l). (8) 

In the equation (8) the derivative ^ is changed with the convectional deriva- 
tive *'^^+(^V) and the term proportional to the pressure gradient VP 
(P is the pressure) is neglected. 

Let us discuss the difference between equation (8) and the equation, usually 
used during the consideration of similar problems in the magnetohydrodyna- 
mical approximation. In previous papers, the acceleration ^ was disregarded. 
The problem was formulated in the inertial frame. It is easy to show, that the 
momentum p in the frame of ZAMOs coincides with the momentum p' in 
inertial frame. Really, according to the definition p^^V , where V—^^, and 

7=«7 (9) 

we will find, that p'=p, because V'—^=^- It is easy then to rewrite the 
equation (8) for the quantities defined in inertial frame. Omitting the prime 
for all quantities, we will obtain 

^+{VV)p=-^aVa+-{E+[VB]). (10) 
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The equation (10) contains the force F=— 7q;Vq!, which was disregarded in 
previous papers. It is the analogy of the force, which causes the centrifugal 
effect and is equal to 



This force exceeds VP {P is the pressure), which was taken into consideration 
in the paper by Kennel et al. (1983). To find the force F, it is not necessary to 
write the equation in the noninertial frame. Staying in the inertial frame and 
taking into consideration that the axially symmetric magnetic field changes 
its direction in the space, i.e. B=B{r)=Br(^-^^ , the equation of motion for 
the guiding center can be obtained. The relativistic generalization of this 
problem (which was stated and solved for the first time by Alfven) can be 
found in Sivukhin (1963). We must mention, that in the paper of Cohen and 
Rosenblum (1972) the problem was investigated in the rotating frame and the 
equations were written for the four-vectors, but then the exact equations for 
the motion of the particles were not used. 

Besides, there are two principal differences of our consideration from pre- 
vious ones (e.g. Michel, 1969; Kennel et al., 1983). In the problem discussed 
here the extremely large surface magnetic field (i?o— 10^^ — 10^^ G) of the neu- 
tron star is essential. On the other hand, in the magnetosphere outside the 
star a magnetic field Bi can be generated. The source of the field Bq is located 
inside the star, but the source of the field Bi is a current j, which is generated 
in the magnetosphere. The current j has no infiuence on the source of the 
pulsar magnetic field Bq and we will consider Bq as an external magnetic field, 
in which the plasma is located. 

Furthermore, in the previous papers another assumption was made, namely 
that the inertia of the particles, i.e. the term which contains the partial time 
derivative ^, may be disregarded. In Henriksen and Norton (1975) the deriva- 
tive ^ for the transversal components of the quantities was changed to the 
derivative with respect to the azimuthal coordinate, and the change 

of the quantities along the magnetic field was disregarded. 

The quantities in equation (10) can be written as: 

E=Eo+Ei, B=Bo+B^, p=Po+Pi, (12) 
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where £"0, Bq and jTq are the basic terms, and Ei, Bi and pi are the pertur- 
bations in the first approximation of expansion over the parameter of weak 
turbulence (-^^). 

^From the condition (1) in zeroth approximation we have: 

Eo+[K)^o]=0 

and from the equation (10) for the radial acceleration we will obtain (Machabeli 
and Rogava, 1994): 



d r Q r 



(13) 



Here we discuss only the radial motion, because the azimuthal one is de- 
termined by corotation (V^ = VLr). The solution of the equation (13) can be 
presented in the form: 

, , Vn SnVlt ,^ 

where Sn and dn are the Jacobian elliptical sine and modulus respectively 
(Abramovitz and Stegun, 1964), Vq is the initial velocity of the particle. Prom 
the equation (13) it follows that if the radial velocity K > the acceleration 
changes its sign and the particle is not accelerated, but braked (see in detail 
Machabeh and Rogava (1994)). 

Using the asymptotic expression for the Jacobian function we find that 
when Vq— >1, 

Vn 

r{t)^-^sinnt. (15) 
For the radial velocity we obtain: 

VQr = VQCOsVtt, (16) 

from which it follows that 

Vo^^{y^,f+{V^^f^const. (17) 

This fact is not quite evident. In the nonrelativistic case, as it is known, 
the particle kinetic energy increases exponentially at the expense of the exter- 
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nal source, which keeps constancy of rotation. This follows from the general 
solution (14) too, when V^-Cl; expanding the Jacobian functions we obtain: 



r{t)=^shnt, (18) 

where sh is the hyperbolic sine. 

The principle difference between nonrelativistic and relativistic motion can 
be easily understood if we rewrite the equation (10) in a form different from 
(13): 

7o^+(7o^^V+7o^)W=7o^^V. (19) 

The second and third terms in the left side describe the change of mass (7- 
factor) in time. It is evident (from (19)) that for small velocities H^— >0 (70— >1) 
we will obtain ^-Q^r=0 and the solution (18). But if (70^00), we 

will obtain ^^+Q^r—0 and the solution (16). Prom the equation (19) also 
follows that the term connected with the change of mass dominates over the 
other terms if 
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Vor>^^^. (20) 

When (20) is valid, the inertia of particles is such, that the centrifugal force 
changes its sign and the motion is braked. We can sec, that when Vq^^I, the 
massless" ("force-free") approximation can be justified only if "force-free" 
means ^=0. If we disregard the right side in the equation (19) in comparison 
with the first term of the left side, this fact automatically means disregarding 
the first term in comparison with the second term in brackets on the left 
side of the equation. Then pQ^=const — the particle motion is uniform and 
straightforward and we have no corotation. Therefore disregarding the force 
F is possible only in the case Q=0. 

Thus the "massless" (" force- free" ) approximation cannot be used for the 
"solid-body" type rotating ('\7=[i1r|) relativistic plasma stream. ^^7^0, espe- 
cially must not be equal to zero. 

If we add to the equation (10) the continuity equation and Maxwell equa- 
tions we shall have a set of equations, which completely describes the behaviour 
of a relativistic plasma stream corotating with the pulsar magnetosphere: 



10 



(VB) = 
(V^) = 47rp 

f = -[V(„^,l 

dE 

— = [V{aB)] - Anaj 

— * 

where p and j are the charge and current densities respectively. 

Conclusion 

We have investigated the role of particle inertia and shown that the "mass- 
less" ("force- free") approximation cannot be used for the "sohd-body" type 
rotating (l/=[Qr]) relativistic plasma stream description. We must mention 
that our investigation was carried out under the assumption £"11 = and if 
£^117^0, this would probably have the strong influence on our results. We also 
derived the equation of motion and a set of equations for the description of a 
relativistic plasma stream in the perpendicular rotator model of pulsar mag- 
netospheres. 
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